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ABSTRACT 


Time dependent variable measures were obtained for all 
major aircraft accidents between July 1971 and July 1974. 
Using these time dependent variables and functional forms 
of these variables, a regression analysis was performed 
for each of eight major aviation commands. By using these 
functional forms of the variables, a relatively high amount 
of variance in aircraft accident rate was accounted for at 
a high confidence level in some eenreanaa When reviewing 
the results of the eight major commands considered, it was 
particularly noted that the variables most instrumental 
in explaining the variance in aircraft accident rate were 
not all pilot oriented but were variables interpreted as 
being related either to pilot experience level, pilot 
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ieee LNTRODUCRRON 


Aviation safety is one of the major tools used to obtain 
a Navy goal of fleet readiness. The reports necessary to 
provide the Naval Safety Center with a means by which to 
judge the Navy's progress in aviation safety are delineated 
in OPNAVINST 3750.6 (Series). Through analysis of reported 
mishaps, the safety center has the ability to determine 
correctable material and personnel deficiencies. Through 
the administrative effort expended in submitting reports, 
the Navy can be repaid many times in the savings of lives 
and material. 

Due to extensive post aircraft accident analysis by 
accident investigation teams, aircraft accidents have been 
broadly categorized in terms of causal factors. Causal 
factors are not to be confused with environmental factors 
but are any event, act, failure/malfunction, circumstance 
or occurrence, the presence or absence of which caused the 
accident. Contributing to accidents are poorly designed 
equipment, supervisory error, improper aircrew training, 
failure of a facility to provide adequate runway drainage, 
maintenance induced equipment malfunctions, etc. 

An accident is designated as a major accident if: 

(1) loss of life occurs; (2) complete loss of an aircraft 


is involved; or (3) substantial damage occurs to any 





aircraft involved where substantial damage is defined in 
Appendix A of OPNAVINST 3750.6 (Series). The most common 
cause cited in the literature on aircraft accidents is 

pilot error. In a study by Brictson, et al. [1969_/, 
approximately seventy-eight percent of the accidents stud- 
ied attributed pilot error as being the primary causal 
factor, with only eight percent being attributed to 

support personnel errors, and the remaining thirteen per- 
cent being attributed to aircraft FaabeS. equipment failure 
and/or weather. 

A study conducted for the Royal Air Force by Goorney 
/ 1965_/, concluded that fatigue, personal worries, com- 
Placency, emotional stress and lack of current flying expe- 
me@ae, directly contributed @o pilot error. Goorney's 
study states that excessive ground duties prior to flying 
attributed to fatigue and that personal worries concerning 
marriages, dating, housing, financial and work problems 
attributed to emotional stress. 

Collicot, et al. / 1972_/ categorized accident causal 
factors into pilot error, material failure, maintenance 
error, and miscellaneous other causes. The authors compared 
Air Force F-4 accidents with Navy-Marine F-4 accidents and 
attributed maintenance error disparities to the fact that 
Air Force F-4 aircraft only realized one-tenth of the can- 
nibalization of parts that Navy-Marine aircraft realized. 
The authors also singled out that the Navy-Marine officer 


job rotational policies resulted in lower in-type flight 





hours for Navy-Marine aviators and thus adversely affected 
their pilot proficiency. The authors also noted that, in 
comparing dual piloted to single piloted aircraft, the dual 
piloted aircraft had fewer accidents per ten thousand hours 
than the single piloted aircraft. The authors conclusion 

was that pilot mental overload was critical in the determina- 
immon of a prime factor in pilot error. 

As a result of the Collicot study and a study by Kowalsky, 
et al. 7 1974_/, the emphasis that pilot induced accidents 
were caused by lack of pilot proficiency was shifted to 
implications that pilot error accidents may be attributed 
to a state of temporary mental overload where the pilots 
incorrectly evaluated information inputed to them during the 
overload period. 

Since the advent of the Naval Safety Center, a great deal 
of effort has been expended in maintaining extensive data 
banks of accident related information. Statistical Analysis 
of available data should enable predictive mathematical 
models to be constructed, as much of the data available can 
be construed as measures of pilot proficiency. Myers / 1974 _/ 
attempted such a statistical analysis, hypothesizing that 
measures of pilot proficiency and experience would suffice 
to form an adequate foundation for accident rate analysis. 

He applied statistical techniques of principle component 
analysis and cluster analysis to ten variables obtained from 
the Individual Flight Activity Reporting System (IFARS) which 


are submitted to the Safety Center. Comparing two groups of 








pilots--one group being involved in aircraft accidents and 
the other group being accident free--his results were not 

as pronounced as desired. The fact that each group con- 
tained only fifty subjects is reason to believe that this 
small sample size suppressed his accident predictive results. 

Maxwell and Stucki / 1975_/ applied regression analysis 
techniques to an all Navy study of aircraft accidents for 
the period July 1968 to June 1974. Their analysis indicated 
that 46.6% of the total accident rate variability is explained 
by pilot related variable measures if one wishes to accept 
a statistical confidence level of 75% as meaningful. Results 
of their study corroborated with the previous studies cited, 
that pilot error of one sort or another is the single largest 
cause of aircraft accidents. 

The author of this writing agrees with Maxwell and Stucki 
that sufficient data should be currently available, from | 
which predictive capability is extractable. The variable 
nature of the monthly accident rate suggests underlying 
factors, causal and thus definable in their role of accident 
perpetration. Using statistical analysis to isolate variable 
measures associated with pilot proficiency and aircraft 
maintenance which vary either directly or indirectly with 
accident rate, then predictive and thus preventive know- 
ledge can assist in lowering the loss of human life and 


the dollar loss resulting from aircraft accidents. 
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If. PROBLEM DEFINITION 


The Naval Safety Center (NSC), Norfolk, Virginia, has 
provided the author with major accident rates computed for 
aircraft accidents for eight major commands: (1) Commander 
Naval Air Forces Atlantic (COMNAVAIRLANT); (2) Fleet Marine 
Forces Atlantic (MARLANT) ; (3) Commander Naval Air Forces 
Pacific (COMNAVAIRPAC); (4) Fleet Marine Forces Pacific 
(MARPAC); (5) Commander Naval Air Training Command (CNATRA); 
(6) Commander Marine Training Command (MARTC); (7) Naval 
Reserves; (8) Commander Research Development Test and 
Evaluation/Naval Air Systems Command (RDT&E/NASC); for 
fiscal years 1969 through 1974. The accident rate is defined 
as the total number of accidents in a given month multiplied 
by a constant factor of-ten thousand and then divided by 
the total monthly hours flown for a given command. Major 
accidents by definition, as delineated in OPNAVINST 3750.6 
(Series), are characterized by loss of life, loss of aircraft, 
or extensive aircraft damage measured in necessary man-hours 
to effect repair. 

A variety of approaches to look for consistent trends 
or cyclic phenomena in aircraft accident rate of Navy-Marine 
accidents was undertaken by Poock / 1976 _/. None were found 
that could pass any statistical tests to verify their 


existence beyond reasonable elements of chance. 


mals 





The purpose of this paper is to further explore accident 
rate dependence on time related variable measures for the 
eight major commands listed above in hopes that one or more 
of these measures can be identified for later use in the 
reduction of accident rates and thus the reduction in the 


loss of lives and money. 
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Itt. ANALYTICAL PROCEDURES 


A. DATA SOURCE 


OPNAVINST 3750.6 (Series) promulgates the requirements 
and procedures for the reporting of each aircraft accident 
or incident involving all Navy and Marine aircraft. This 
establishes control over accident data with the goal of 
increasing aviation safety. The reports are forwarded 
to the Naval Safety Center (NSC) for inclusion in their 
master data bank. Accident data currently available from 
NSC can provide approximately eighty variable measures for 


each accident. 
B. DATA SELECTION 


The NSC data bank provides a ready source of data for 
each accident or data set. However, as this is a continuation 
of the work of Maxwell and Stucki, the 2110 computer data 
cards obtained from NSC for their study were also used in 
mits study. 

The initial step in the conduct of the current accident 
rate analysis was to select appropriate variable measures or 
data points. A data point for an accident was considered 
to be any suitable variable measure associated with the 


accident. and a data set consisted of data points for a 


1, 





specific accident. Specific measures were chosen in mutual 
discussions with NSC personnel. 

selection of appropriate data points required that each 
point be time dependent. Data point time dependency and 
Subsequent selection was based on the variable descriptions 
contained in the Manual of Code Classification for Navy 
Aircraft Accident, Incident and Ground Accident Reporting 
(Code Manual) promulgated by NSC. 

A sufficient number of data sets had to be incorporated 
into the analysis to facilitate viable statistical results, 
but the span of time defined by the data sets had to be 
chosen with care. Unfortunately, of the six fiscal years 
of data available to the author, only fiscal years 1972, 
1973, and 1974 were used in this study. Fiscal years 1969, 
1970, and 1971 were eliminated primarily because some data 
sets were either incomplete or entirely void of information. 
A change in reporting procedures during 1972 alleviated this 
problem. Also, in the fiscal years chosen for the analysis, 
a particular accident too often had to be disregarded due to 
incomplete data or to the fact that the particular occurrence 
involved another aircraft and the other aircraft was deemed 
at fault. Therefore, either due to insufficient data or 
to an occurrence and not an accident, a total of 636 of 
the original 2110 data sets available were used. Data for 
each variable of these 636 data sets were complete and there 
is absolutely no data missing in the data sets used in 


the analysis. 
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Ten data points from each set of the original thirty 
data points requested by Maxwell and Stucki were selected 
for inclusion in the current analysis. Nine of these data 
points are the same data points used in the Maxwell and 
Stucki study with the addition of the data point "Years 
Experience as a Designated Naval Aviation" and the data 
point "Flight Purpose Code" deleted. In addition, the data 
points selected also allowed the author to construct four 
additional variables deemed as measures of experience and 
pilot proficiency. Table I lists the basic variables used 
in the study. The first ten variables listed are from the 
Naval Safety Center and variables eleven through fourteen 
are the constructed variables. 

Pilot age, years experience as a designated Naval aviator, 
experience, age at designation as a Naval aviator and total 
flight time in the aircraft involved in the reported accident 
have been considered as measures of pilot experience. If 
these variables are true measures of experience they should 
exhibit negative correlation with accident rate, assuming 
a relationship exists. 

There is the age old adage of "practice makes perfect" 
which this author construes as a measure Of (pilot preumeleney. 
Included in the category of pilot proficiency were total 
flight time during the preceding ninety days, total night 
flight time during the preceding ninety nights, total day 
flight time during the preceding ninety days, the number 
of carrier landings in the last thirty days, the number of 


carrier landings in the last thirty nights, and the ratio 
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TABLE I 
DATA SET INCLUDED IN CURRENT STUDY 
Accident rate by month (RATE) 
Pilot's age (AGE) 
Years experience as a designated Naval aviator (DNA) 


Total flight time in accident involved aircraft model 
(TTIME) 


Total flight time during preceding ninety days (TOT90) 


Total night flight time during preceding ninety nights 
(NITE9O) 


Daylight carrier landings during preceding thirty days 
(CLDAY ) 


Night carrier landings during preceding thirty nights 
Goran it) 


Number of aircraft tours (ACTOUR) 


Aircraft flight hours since last major or minor 
inspection (ACHRS) 


The ratio of years experience as a designated Naval 
aviator and total flight time in the accident involved 
aircraft model (EXPER=DNA/TTIME ) | 


Age at designation as a Naval aviator (WINGS=AGE-DNA) 
Total flight time during preceding ninety days minus 
total night flight time during preceding ninety nights 
(DAY90=TOT90-NITE9O ) 

The ratio of total night flight time during the pre- 


ceding ninety nights to total day flight time during 
the preceding ninety days (NITEDAY=NITE90/DAY90) 
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of night flight time to day flight time during the preceding 
ninety days and nights. 

To measure airframe age and general condition, two var- 
iables were selected. All Navy/Marine aircraft are required 
to undergo a Periodic Aircraft Rework (PAR) cycle for analysis 
and repair after having accumulated a specific number of 
flight hours. Thus, the number of aircraft tours was chosen 
as one variable. As a measure of aircraft condition, the 
aircraft flight hours since the last major or minor inspection 


was also selected as a variable. 
C. PRELIMINARY DATA PREPARATION 


A common theoretical proposition in parametric statistics 
states that changes in one variable can be explained by 
reference to changes in several other variables. Such a 
relationship is described in a simple way by a multiple 
linear regression equation as described in Appendix C. To 
use this statistical method the relationships between the 
variables require the following assumptions: (1) normality; 
(2) data must be in the interval measurement scale; 

(3) the number of observations exceeds the number of coef- 
ficients to be estimated; and (4) no exact linear relation- 
ships exist between any of the explanatory variables. 
Normality was assumed by invoking the Central Limit Theorem ° 
and the technique of averaging by month was used to transform 


the data into interval data. 


Ie 





Raw data for each of the first ten variables of Table I 
was averaged by month for each of the thirty-six months within 
the time span selected and was also used to compute the values 


for the remaining variables of Table I. 
D. THE ANALYSIS TECHNIQUE 


To conduct the statistical analysis of the data set the - 
author selected the forward (stepwise) multiple regression 
computer program package developed by Jae-On Kim and Frank 
J. Kahout at the University of Iowa. The program is included 
in the Statistical Package for the Social Sciences (SPSS) 
compiled and edited by Nie, et al. / 1975_/. 

Multiple regression may be viewed as a descriptive tool 
whereby the linear dependence of one variable on other var- 
lables is summarized and decomposed. Kim and Kahout state 
that the most important uses of the technique as a descriptive 
tool are: 

(1) to find the best linear prediction equation 

and evaluate its prediction accuracy; (2) to 

Comino. for.ouner contoundins factors in order 

to evaluate the contribution of a specific 

variable or set of variables; and (3) to find 

structural relations and provide explanations 

for seemingly complex multivariate relation- 

ships, such as is done in path analysis. 

However, the primary purpose is to evaluate and measure over- 
all dependence of a specified (dependent) variable on a set 
of other variables (independent). The dependent variable 
used for the current study was accident rate and the inde- 


pendent variables consisted of those variables listed as 2 


through 14 in Table I and all the variables listed in Table II. 
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To explain the variables listed in Table II it must 
again be stated that regression analysis oo the under- 
lying relationships among the variables are linear and addi- 
tive. There are many occasions for which such simple linear 
models are inadequate. Early efforts by the author dealt 
entirely with first order variables but examination of 
scatterplots of the Peace which are conceived as meas- 
ures of the error component, indicated lack of linearity and 
the presence of curvilinearity. Therefore, functional terms 
were added to the variable list as suggested by Kim and Kahout. 
These additional variables are listed in Table II. 

The SPSS computer program is designed to provide the user 
with a considerable number of control options. The listwise 
deletion option was chosen for use by the author as it is 
the most conservative and accurate of the OpELONs. 

The forward (stepwise) multiple regression technique 
(Appendix C) is particularly useful in studies of the current 
type. It is appropriate to enter independent variables one 
by one on the basis of a pre-determined statistical criteria. 
This procedure is used when a researcher desires to isolate 
a subset of available predictor variables that will yield an 
optimal prediction equation with as few terms as possible. 

Draper and Smith / 1966 _/ conclude the following in their 
discussion of various regression procedures: "We believe 
this (stepwise regression) to be the best of the variable 
selection procedures discussed and recommend its use." 


The percentage of the accident rate variability explained 


ey 





Oo On WAN FWD FH 


PPP PPP PRP eR 
0 ON AWN FWD KB 


TABLE II 


CONSTRUCTED VARIABLES 


AGE2 = (AGE) © 

TTIME2 = (TTIME)* 
TOT902 = (TOT90)“ 
WINGS2 = (WINGS) © 
DAY902 = (DAY90)“ 
CLDAY2 = (CLDAY)“ 
CLNIT2 = (CLNIT)¢ 
ACHRS2 = (ACHRS)“ 
EXPER2 = (EXPER) * 
DNA2 = (DNA)“ 


NITE902 = (NITE90)“ 
ACTOUR2 = (ACTOUR) 
NITEDAY2 = (NITEDAY) “ 
PETEDAYT = 1.0/NITEDAY 2 
Peet — 0 /AGE2 

Mire = 1.6/TTIME2 
TOTQOOI = 1.0/TOT902 
NITE9OI = 1.0/NITE902 
DAawGOl =n1\,.0/DAY9O2 


Z0n 
2A 
Le 
Dae 
24, 
25. 
26. 
oe 
Don 
29. 
30. 
31. 
32. 
Sit 
34. 
J 
30. 
37. 
38. 
Bos 


20 


CLEAYT = 0, cupay2 
Chine —"4. 0/7 Chir Ts 
ACTOURI = 1.0/ACTOUR2 
ACHRSI = 1.0/ACHRS2 
DNAI = 1.0/DNA2 
EXPERI = 1.0/EXPER2 
WINGSI = 1.0/WINGS2 
RTAGE = (AGE) * 

RTDNA = (DNA)? ; 
RTTTIME = ( TTIME) * 
REGO Cm (TOT9O) * 
RTDAY9O = (DAY90) * 
Reno — (CLDAY) * 
RICINIT = (CLNIT) * 
RTACHRS = (ACHRS) * 
RTEXPER = (EXPER) * 
RTWINGS = (WINGS)* ; 
RINITE9O = (NITE9O) * 
RTACTOUR = (ACTOUR)? : 
RTNITEDAY = (NITEDAY)? 





by the time related independent variables entered is exactly 
the type of statistical output required to ascertain the 
causal factors responsible for that variability. 

Thirteen separate regression calculations were attempted 
for each of the eight commands considered. Each of the 
separate regression calculations were permutations of the 
variables listed in Table I and the variables listed in 
Table II. Regression calculations using only the first 
eight variables of Table I or regression calculations using 
the first ten variables of Table I led to the suspicion of 
curvilinearity in the data as a result of an inspection of the 
residual scatter plots. Therefore, the author attempted 
various regression calculations using combinations of the 
variables of Table I with either the squares of those var- 
iables, the inverse squares of those variables, the square 
roots of those variables or any combinations thereof. 

There is no unique procedure for selecting the best 
criteria and personal judgement is a necessary part of any 
decision reached by the analyst. The author set up guide- 
lines to choose the "best" equation as follows: (1) the 
regression equation would be limited to five or six inde- 
pendent variables; (2) the regression equation chosen should 
account for the highest amount of the variance in aircraft 
accident rate; (3) the residual plots would be examined 
for the best indication of minimum residuals; and (4) the 
statistical significance of each regression would be cal- 


culated as outlined in Appendix D and the regression equation 


eu 





yielding the greatest amount of accountable variance with 


the highest confidence level would be the equation selected. 
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IV. RESULTS 
A. COMNAVAITRLANT 


Early efforts at regression analysis by the author used 
only the variables listed in Table I. These efforts led to 
ae SUSPIC1ON Of curvilinearity as a result of examination of 
the residual plots. The order of variable inclusion at the 
early stage was: (1) total flight time during the ninety 
days preceding the alec uae OmoT90); and (2) aimerani er meen, 
hours since the last major or minor inspection (ACHRS). In 
equation form the regression became: 


Bate = 0.952905 + On00653)°70T90) +°0.00163( ACHES) 


This equation accounted for 19.5% of the variance Phy a) ieee pecan 
accident rate at a 95% confidence level. 

Final regression results using functional forms of the 
variables indicated that the hierarchical order of variable 
inclusion as governed by the individual variable contributions 
towards explaining accident rate variance was: (1) the square 
root of aircraft flight hours since the last major or minor 
inspection (RTACHRS); (2) the square of aircraft flight 
hours since the last major or minor inspection (ACHRS2) ; 

(3) the inverse square of the total flight time during the 
ninety days preceding the accident (TOT90I); (4) the inverse 
square of the quotient years experience as a designated 


Naval aviator and total flight time in the accident involved 
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aireraft model (EXPERT); and (5) thessequearce roo. of tac 
quotient years experience as a designated Naval aviator 
and total flight time in the accident involved aircraft 
model (RTEXPER). Table III lists the simple correlation 
coefficients and Table IV is a summary listing of computer 
output provided by the SPSS package. In equation form the 


regression became: 


RATE = 0.89771 + 0.11823(RTACHRS) - 0.00002(ACHRS2) 
- 464,21653(TOT90I) - 0.00002( EXPERT ) 
- 4,76996( RTEXPER) 


or in terms of Table I variables 


i 
2 


RATE = 0.89771 + 0.11823(ACHRS)2 - 0.00002(ACHRS) © 


-~ 464, 21653(TOT9I0)~ - 0.00002(TTIME/DNA) < 


1 
2 


-~ 4,76996(DNA/TTIME) 


These equations account for 42.9% of the variance in aircraft 
accident rate at a 99% confidence level. Therefore, by 
using functional forms of the variables, an additional 
23.4% of the variance was accounted for--a substantial 
improvement in the results. 

The net effect of the functional forms of hours since 
the last major or minor inspection is positive. The remain- 
ing variables are interpreted as being measures of either 
pilot experience or pilot proficiency and they have a net 


effect of reducing aircraft accident rate. 
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B. MARLANT 


Using the same early methods as those used for COMNAVAIRLANT 
the order of variable inclusion at the early stage was: 
(1) years experience as a designated Naval aviator (DNA); 
(2) daylight carrier landings during the preceding thirty 
days (CLDAY); (3) the quotient of years experience as a 
designated Naval aviator and total flight time in the accident 
involved aircraft model (EXPER); (4) the quotient of night 
flight hours to day flight hours during the preceding ninety 
days (NITEDAY); and (5) total flight time in accident 
involved aircraft model. In equation form the regression 


became: 


RATE = 0.96753 + 0.09865(DNA) + 0.08248(CLDAY) 
- 20.27987(EXPER) + 1.41098(NITEDAY ) 
- 0.00048(TTIME) 
This equation accounted for 32.6% of the variance in air- 
craft accident rate at a 75% confidence level. 

Final regression results using functions of the var- 
jables indicated that the hierarchical order of variable 
inclusion as governed by the individual variable contri- 
butions towards explaining accident rate variance was: 
(1) the square root of years experience as a designated 
Naval aviator (RTDNA); (2) the square root of total 
night flight time during the preceding ninety nights 
(RINITE9O); (3) total night flight time during the 
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preceding ninety nights (NITE90); (4) the quotient of 
years experience as a designated Naval aviator and total 
flight time in the accident involved aircraft model (EXPER); 
fos) total flightietime in the accident mm .o h.edmomer-ceum 
model (TTIME); and (6) the square of night carrier land- 
ings during the preceding thirty nights (CINIT2). Table V 
lists the simple correlation coefficients and Table VI is a 
summary listing of computer output provided by the SPSS 


package. In equation form the regression became: 


RATE = -0.17128 + 0.88657(RTDNA) + 0.52634(RTNITE9O) 
-0.06351(NITE90) - 57.37027(EXPER) - 0.00115(TTIME) 
SGU eGilul Valetta) 


Or in terms of Table I variables 


i 
2 


RATE = -0.17128 + 0.88657(DNA)® + 0.52634(NITE9O) 
- 0.06351(NITE9O) - 57.37027(DNA/TTIME) 
~ 0.00115(TTIME) + 0.21877(CLNIT)“ 


These equations account for 56.5% of the variance in air- 
craft accident rate at a 99% confidence level versus the 
32.6% accounted for initially. A substantial improvement 
was again observed when using functional forms of the 
variables in the regression analysis. 

When reviewing the correlation matrix (Table V) it is 
noted that every variable present in the equation is posi- 
tively correlated with accident rate except for EXPER (a 


measure of pilot experience) and the square of night 


carrier landings in the last thirty days (a measure of pilot 
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proficiency). The net affect of the functional forms of 
night flight hours during the last ninety nights is positive, 
an indication that night flying may be a hazardous evolution 


and add to accident rate. 
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C. COMNAVAIRPAC 


Initial efforts yielded the following order of variable 
inclusion: (1) total flight time during the preceding 
ninety days (DAY90); (2) the number of aircraft tours 
(ACTOUR); and (3) years experience as a designated Naval 


aviator. In equation form the regression became: 


RATE = -0.07944 + 0.01289(DAY90) + 0.10561(ACTOUR) 
+ 0.04137(DNA) 
This equation accounted for 33.8% of the variance in air- 
craft accident rate at a 99% confidence level. 

Final regression results using functions of the var- 
lables indicated that the hierarchical order of variable 
inclusion as governed by the individual variable contri- 
butions towards explaining accident rate variance was: 
(1) the inverse square of total day flight time during 
the preceding ninety days (DAY90I); (2) the square of 
the age of the pilot at designation as a Naval aviator 
(WINGS2);: (3) the number of aircraft tours (ACTOUR); 
(4) the inverse square of aircraft hours since the last 
major or minor inspection (ACHRSI); and (5) the square 
of the hours since the last major or minor inspection 
(ACHRS2). Table VII lists the simple correlation coef- 
ficients and Table VIII is a summary listing of computer 
output provided by the SPSS package. In equation form the 


regression became: 
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RATE = 1.47479 - 670.77272(DAY90I) - 0.00134(WINGS2) 
+ 0.15733(ACTOUR) + 304.94442(ACHRST) 
+ 0.00001(ACHRS2) 


wenn terms of Table 1 variables 
RATE = 1.47479 - 670.77272/(DAY90)* - 0.00134(AGE - DNA) 


+ 0.15733(ACTOUR) + 304.94442(AcHRS) - 


+ 0.00001(ACHRS) “ 


These equations account for 45.7% of the variance in air- 
craft accident rate at a 99.5% confidence level, an increase 
of 11.9% in accounting for variance. 

It is noted in these equations that the number of 
aircraft tours and functional forms of hours since the 
last major or minor inspection have a positive effect on 
accident rate. Conversely, variables interpreted as 
Medeurés of pilot proficiency and experience level have 


a negative effect. 
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D. MARPAC 


Initial efforts at regression analysis yielded very poor 
results for this command. The only variable to be included 
in the equation was the quotient of night flight hours to 
day flight hours during the preceding ninety days. In 


equation form the regression became: 
RATE = 1.75337 - 2.47720(NITE90/DAY90) 


This equation accounted for 9.2% of the variance in air- 
Sratt accident rate at a 75% confidence level. 

Further regression efforts using functions of the var- 
lables led to much better results with the hierarchical 
Bader Of variable inclusion being: (1) the square of 
the number of aircraft tours (ACTOUR2); (2) the number 
of aircraft tours (ACTOUR); (3) the square of total 
night flight time during the preceding ninety nights 
(NITE902); (4) the square of the quotient total night 
meet) T nours to total day flight hours during the preceding 
Patety days (NITEDAY2): (5) aircraft flight hours since 
the last major or minor inspection (ACHRS); and (6) the 
Square of age at designation as a Naval aviator (WINGS2). 
Table IX lists the simple correlation coefficients and 
Table X is a summary listing of computer output provided 
by the SPSS package. In equation form the regression 


became: 
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RATE = 1.26335 - 0.32047(ACTOUR2) + 2.21178(ACTOUR) 
+ 0.01042(NITE902) - 19.20769(NITEDAY2) 
- 0.00565(ACHRS) - 0.00462(WINGS2) 


or in terms of Table I variables 


RATE = 1.26335 - 0.32047(ACTOUR)* + 2.21178(ACTOUR) 
+ 0.01042(NITE90)“ - 19.20769(NITE90 /DAY90) @ 
- 0.00565(ACHRS) - 0.00462(AGE-DNA)* 


Either equation accounts for 39.2% of the variance in air- 
craft accident rate at a 95% confidence level, a substantial 
improvement over the 9.2% of variance accounted for initially. 
Reviewing the correlation matrix it is noted that every 
variable present in the equation is negatively correlated 
with accident rate. However, the intercorrelations of the 
variables has the affect of causing some of the signs of 
the coefficients in the equation to change, with the net 
eect of the functional forms of alrcrafit tours being 


positive. 
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E. CNATRA 


All attempts at regression analysis for this command 
yielded essentially the same results. Early efforts using 
only the variables listed in Table I yielded the following 
order of variable inclusion: (1) the quotient of total 
moms Llirght hours so total day flight hours durime the 
preceding ninety days (NITEDAY); (2) total flight time 
in accident involved aircraft model (TTIME); (3) the 
mmertent of total flight time in.accident involved alrcrarft 
model to years experience as a designated Naval aviator 
(TTIME/DNA); (4) aircraft flight hours since the last 
major or minor inspection (ACHRS); and (5) the quotient 
of years experience as a designated Naval aviator to total 
flight time in accident involved aircraft model. In equation 


form the regression became: 


RATE = 5.22976 + 51.98496(NITE90/DAY90) + 0.00955( TTIME) 
- 0.01973(TTIME/DNA) - 0.04198(ACHRS) 
- 216.41848(DNA/TTIME) 


imeneauation accounted) fom 32.3% of the variance in air= 
craft accident rate at a 90% confidence level. 

Later regression results using functions of the variables 
yielded: (1) the square root of aircraft flight hours since 
the last major or minor accident (RTACHRS); (2) the square 


of aircraft hours since the last major or minor inspection 
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(ACHRS); (3) the square of the total day flight hours 
flown during the preceding ninety days (DAY902); (4) the 
Square root of the quotient of total flight time in the 
accident involved aircraft model to years experience as a 
designated Naval aviator (1.0/RTEXPER); and (5) the 
quotient of total flight time in accident involved air- 
craft model and years experience as a designated Naval 
aviator. Table XI lists the simple correlation coefficients 
and Table XII is a summary listing of computer output pro- 
vided by the SPSS package. In equation form the regression 


became: 


RATE = -0.31346 + 0.11248(RTACHRS) - 0.00003(ACHRS2) 
- 0.00002(DAY902) + 0.03482(1.0/RTEXPER) 
~ 0.00092(1.0/EXPER) 


or in terms of Table I variables 


a 
RATE = - 0.31346 + 0.11248(ACHRS)2 - 0.00003(ACHRS) “ 
2 
~ 0.00002(DAY90)* + 0.03482(TTIME/DNA) 2 
- 0.00092( TTIME/DNA) 


Either equation accounts for 33.0% of the variance in air- 
craft rate at a confidence level of 90%; representing no 
significant change in the results. The functional forms 
of aircraft hours since the last major or minor inspection 
and the reciprocal forms of EXPER have a net positive 
effect while the square of day flight time during the last 


ninety days has a negative affect. 
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FP. NAVAL RESERVES 


Initial efforts yielded the following order of variable 
inclusion: (1) total night flight time during the preceding 
ninety nights (NITE90); (2) the quotient of total night 
flight time and total day flight time during the preceding 
ninety days and nights (NITE90/DAY90); (3) total flight 
time during the ninety days preceding the accident (TOT90); 
and (4) the number of daylight carrier landings during the 
preceding thirty days. In equation form the regression 


became: 
RATE = 0.97395 + 0.06986(NITE9O) - 1.22954(NITE9O/DAY90) 
= ©.00777( TOTSO) - GRO3323(CLDAY) 


This equation accounted for 71.4% of the variance in air- 
craft accident rate at a 99% confidence level. 

Final regression results using functions of the variables 
indicated that the hierarchical order of variable inclusion 
as governed by the individual variable contributions towards 
explaining aircraft accident rate variance was: (1) the 
square of total night flight time during the preceding 
ninety nights (NITE902); (2) the square of the quotient 
of total night flight time and total day flight time during 
the preceding ninety days (NITEDAY2); (3) the inverse 
Square of aircraft flight hours since the last major or 
minor inspection (ACHRSI); (4) the square root of daylight 


carrier. landings during the preceding thirty days; and 
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(5) the inverse square of the quotient of years experience 

as a designated Naval aviator and total flight time in the 
accident involved aircraft model (EXPERI). Table XIII lists 
the simple correlation coefficients and Table XIV is a summary 
listing of computer output provided by the SPSS package. In 


equation form the regression became: 
RATE = 0.83553 + 0.00084(NITE902) - 0.56367(NITEDAY 2) 


-311.23233(ACHRSI) - 0.22889(RTCLDAY) 
+ 0.00002(EXPERT) 


Or in terms of Table I variables 


mee =) 0.383553 + 0.00084(NITE9O) “ = 0. 56367(NITE90/DAY90) “ 


1 
2 


~ 311.23233(1.0/ACHRS)* - 0.22889(CLDAY) 
+ 0.00002(TTIME/DNA) © 


Either equation accounts for 81.3% of the variance in air- 
craft accident rate at a 99.9% confidence level, a 9.9% 
increase of accountable variance over initial efforts. The 
net effect of the square of night flight hours during the 
preceding ninety nights and the quotient of night flight 
hours and day flight hours during the preceding ninety 
nights and days is negative. This may indicate that the 
meehner the ratiwo of ment flight hours: teeday pement heuns 
(a measure of pilot proficiency), the lower the accident 


rate. 
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G. MARTC 


Tnetial efforts using only thesvariablecmor = ables 
yielded the following order of variable anentueucnty (1) the 
Mulocrenb Of “Gotal night flight hours and Gotal day firens 
hours during the preceding ninety days (NITEDAY); (2) total 
flight time in accident involved aircraft model (TTIME) ; 
Meeeune quotient of total flight time in accidengeinvorvcad 
aircraft model and years designated as a Naval aviator 
(TTIME/DNA); (4) age at designation as a Naval aviator 
Mmnes): and (5) number of aircraft tours (ACTOUR). In 


equation form the regression became: 
RATE = - 45,84436 + 32.23927(NITEDAY) + 0.02707(TTIME) 


~ 0.15004(1.0/EXPER) + 2.20096(WINGS) 
~ 0.75295(ACTOUR) 


This equation accounted for 68.0% of the variance in air- 
craft accident rate at a 95% confidence level. 

Final regression results using functions of the variables 
indicated that the hierarchical order of variable inclusion 
as governed by the individual variable contributions towards 
explaining aircraft accident rate varlance was: (1) the 
Square of the quotient total night flight time and total 
day flight time during the preceding ninety days (NITEDAY 2) ; 
(2) total flight time in accident involved aircraft model 
(TTIME); (3) the inverse square of the quotient years expe- 


rience as a designated Naval aviator and total flight time in 
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the accident involved aircraft model (EXPERI); (4) total 
night flight time during the preceding ninety nights (NITE9O) ; 
and (5) the square of age at designation as a Naval aviator 
(WINGS2). Table XV lists the simple correlation coefficients 
and Table XVI is a summary listing of computer output provided 


by the SPSS package. In equation form the regression became: 


RATE = - 0.36132 + 65.28325(NITEDAY2) + 0.00328(TTIME) 
- 0.00008(EXPERI) - 0.41913(NITE9O) 
+ 0.00385(WINGS2) 


or in terms of Table I variables 


RATE = - 0.36132 + 65.28325(NITE90/DAY90)* + 0.00328(TTIME) 
- 0.00008(TTIME/DNA)* - 0.41913(NITE9O) 
+ 0.00385(AGE-DNA) © 


Either equation accounts for 79.9% of the variance in air- 
craft accident rate at a 95% confidence level, a 10.9% increase 
of accountable variance over initial efforts. 

It should be noted for this command that only one variable 
interpreted as a measure of experience (EXPERI) is negatively 
correlated with accident rate. Another measure of experience 
(WINGS2) and three variables interpreted as measures of pilot 
proficiency are all positively correlated with accident rate. 
This positive correlation of pilot proficiency and experience 
variables is contrary to what the author hypothesized. The 
Study by Poock / Roo. 7, also concluded that this command 


produced unusual results. 
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H. RDT&E/NASC 


These two commands were grouped together due to the fact 
that neither command had enough data points to be analyzed 
if considered as individual commands. Initial efforts using 
only the variables of Table I yielded the following order of 
marteaole inclusion: (1) total night flight hours cee the 
preceding ninety nights (NITE90); (2) the quotient of 
years experience as a designated Naval aviator and total 
flight time in accident involved aircraft model (EXPER) ; 
and (3) pilot's age. In equation form the regression 


became: 


RATE = 0.93622 - 0.08430(NITE90) + 0.31653(EXPER) 
- 0.01623( AGE) 


This equation accounted for 56.5% of the variance in air- 
craft accident rate at a 95% confidence level. 

Final regression results using functions of the variables 
indicated that the hierarchical order of variable inclusion 
as governed by the individual variable contributions towards 
explaining aircraft accident rate variance was: (1) the 
square of total night flight during the preceding ninety 
nights (NITE902); (2) the square root of total flight 
'time in accident involved aircraft model (RTTTIME); (3) the 
inverse square of years experience as a designated Naval 


aviator (DNAI); (4) the square root of the inverse quotient 


- 





of years experience as a designated Naval aviator and 

total flight time in the accident involved aircraft model 
(1.0/RTEXPER); (5) the inverse square of total flight 

time during the preceding ninety days (TOT90I); and (6) the 
Square root of total flight time during the preceding ninety 
days (RTTOT90). Table XVII lists the simple correlation 
coefficients and Table XVIII is a summary listing of com- 
puter output provided by the SPSS package. In equation 


form the regression became: 


RATE = - 1.21331 + 0.71262(NITE902) + 0.23598(RTTTIME) 
+ 19.68521(DNAI) - 0.61421(1.0/RTEXPER) 
+ 182.91525(TOT9OL) + 0.38527(RTTOTIO) 


or in terms of Table I variables 


Lt 
a 


RATE = - 1.21331 + 0.71262(NITE90)* + 0.23598(TTIME) 


= 
2 


+ 19.68521(1.0/DNA) - 0.61421(TTIME/DNA) 


a 
= 


+ 182.91525(1.0/TOT90)* + 0.38527(TOT9O) 


Fither equation accounts for 80.8% of the variance in air- 
craft accident rate at a 95% confidence level, a 25.3% 


mmcrease Over initial efiorts. 
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V. DISCUSSION 


Reviewing the results it can be seen that as functional 
forms of the variables were introduced into the regression 
analysis that more variance in aircraft accident rate was 
accounted for at a higher confidence level. This was true 
for all eight commands considered. 

Observing the signs of the correlation coefficients 
listed in Table XIX, it can be seen that the variable EXPER, 
which is the ratio of years experience as a designated 
Naval aviator and total flight time in accident involved 
aircraft model, was negatively correlated with accident 
rate in all cases. This indicates that the more flight 
time in a particular aircraft model that a pilot possesses 
per year of designated Naval aviator service that he will 
become more safe and thus have fewer accidents. This could 
suggest that the more experienced a pilot becomes the larger 
a repertoire of near tragedies he has to draw from and thus 
the more reminders or analogies he has to compare with 
current situations. 

Conversely, the variable CLNIT, which is the number of 
night carrier landings a pilot has during the last thirty 
nights, showed positive correlation with aircraft accident 
rate in the three commands where the variable was present 


(COMNAVAIRLANT, MARLANT, COMNAVAIRPAC). This may indicate 


48 





1 + + + 


+ + + 


+ 


ee 


OSVN 
/ ASL 


+ + + + 


Ob VIN 


+ + + 1 


+ + + 


SHUAVN 


- 
+ + + 1 
+ ++ + + + + 
+ + + + + + 


I 

a 

+ + 
+ + + 


VabVNoO OVdYd VN OVdY IV DNVIYVN DNYIYLV 


AVNINOD 


SUNATOLHHHOO NOLILVIEYYOO AO SNDIS 


XIX WIdVb 


AVNINOO 


AVGULIN 
SdHOV 
qnobov 
GLINTIO 
AVA'TO 
O6AV0 
O6aLIN 
O6L0L 
ddd Xa 
SONIM 
HIND Dib 
VNd 
GOV 


aTaVIYVA 


49 








that night carrier landings are more taxing evolutions 
andevend Gor add@ tema reratiaadeecuden perace. 

It has been common belief that pilot proficiency lessens 
aircraft accident rate. However, when again reviewing the 
correlation coefficients, it is observed that the correlation 
of total flight time during the preceding ninety days (TOT90) 
and day flight time during the preceding ninety days (DAY90) 
is positively correlated with aircraft accident rate except 
in the cases of CNATRA and MARTC (both training commands). 

It is also observed that the variable night flight time 
during the preceding ninety nights is positively correlated 
with aircraft accident rate in all commands except CNATRA. 
Considering that the training commands are on a schedule 
where the pilot should not experience fatigue, the positive 
correlation of these variables with aircraft accident rate 
would tend to indicate that if flight time is not uniformly 
spread out that practice does not necessarily make perfect 
and that a possible fatigue factor may be causing an increase 
in aircraft accident rate. Fatigue could be construed as 
having the majority of flight time during the preceding 
ninety days massed into a short period of time instead of 
being distributed uniformly over the ninety day period. (It 
is understood that data for a thirty day period is now being 
recorded.) Another explanation would support Goorney's 
supposition that pilot complacency may increase directly 

as the number of hours flown and thus contribute to aircraft 


accidents. In this case, new pilots in the training command 


DO 





most likely will not become complacent or they will not 
receive their wings. 

The remaining variables listed in Table XIX did not 
indicate a dominant sign of the correlation coefficient 
common to all commands. Therefore, the author will not 
attempt to explain the remaining variables individually. 

When reviewing the results of the major commands con- 
sidered, it should be noted that the variables most instru- 
mental in explaining the variance in aircraft accident 
rate are not all pilot oriented variables as was shown in 
the Maxwell and Stucki study which was an overall Navy/ 
Marine study. The overall results don't seem to apply 
when taking a more microscopic look at individual commands. 
COMNAVAIRLANT, COMNAVATRPAC, MARPAC, CNATRA and Naval Reserves 
regression equations all contain variables interpreted as 
being either related to experience level, pilot proficiency, 
or aircraft condition. Only MARLANT, MARTC, and RDT&E/NASC 
meceunt tor the varlance Im alrerart accident rate by 
variables interpreted as being pilot oriented variables 
in their respective regression equations. Condition of 
the aircraft must therefore also be considered when attempt- 


ing to account for the variance in aircraft accident rate. 


oe 





VI. RECOMMENDATIONS 


The current study includes eight major commands in its 
treatment of aircraft accident rate. Functional forms of 
the variables produced the best results. For this type of 
study however, the author recommends that the distributions 
for each variable be ascertained and that these distributions 
whether they be gamma, beta, exponential, etc. be used in 
an attempt to analyze aircraft accident rate. 

This study was limited to consideration of only accident 
mivolved pilots and accident involved aircraft. A future 
study could be done comparing these pilots aa aincrany 
to accident free pilots and aircraft in order to ascertain 
whether these groups are from the same distribution. It 
would then be possible to ascertain critical points ina 
erlot’s career or the life of an aircraft. 

A more microscopic study than this study may seem appro- 
priate. However, subdividing the data any further into 


Airwings may limit the analyst to only a small number of 


data points. 
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APPENDIX C 
FORWARD (STEPWISE) MULTIPLE REGRESSION 


The basic principles of regression analysis may be 
extended to situations involving two or more independent 
variables. The general mathematical form of the (unstand- 


ardized) regression is 


Y'=Art Sl a BAX a canes Be Xe 


where Y* represents the estimated value of the dependent 
variable, A is the Y intercept, the B, are regression 
coefficients, and the X; are the ee eence nc variables. 
It is assumed that this is a complete set of variables 
from which the equation is to be chosen and includes any 
functions such as squares and inverse squares thought to 
be desirable and necessary. 

The A and B; coefficients are selected in such a way as 
to minimize the sum of squared residuals, d(x-xy?. By 
minimizing the squared residuals, the regression technique 
maximizes the correlation between the actual dependent 
variable (Y) and the estimated dependent variable (Y') 
while the correlation between the independent variables 
and the residual values (Y-Y') are reduced. 

The proportion of variance of Y explained, i.e., the 
goodness of fit of the regression equation can be evaluated 
by examining the square of the multiple correlation Qa 


where Ro is calculated by: 
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Se 
— t —— 
2 55 59 a2 DD zs DY i 


DD 59 (Y i) (Y nt) 
Or 


variation in Y explained by the combined 
“5 linear influence of the independent variables 





moval Variation an Y 


where Ss, is the total variation or sum of squares in Y, 
aa is the sum of squared residuals, and SS neg is the 
regression of squares. 

The forward (stepwise) selection procedure available in 
the Statistical Package for the Social Sciences (SPSS ver- 
sion 6) used for this study inserts variables in turn until 
the regression equation is satisfactory. Independent var- 
lables are entered only if they meet certain statistical 
Gmimeria, The order of inclusion is determined by the 
respective contribution of each variable that explains the 
greatest amount of variance previously unexplained by the 
variables already in the equation. The variable that does 
this is the next variable to be entered. 

This order of insertion is determined by using the 
partial correlation coefficient as a measure of the impor- 
tance of variables not yet in the equation. The, basic 
procedure is as follows. First select the X most correlated 
with Y (suppose it is X,) and find the linear regression 


A 
equation Y = £(X,). Next find the partial correlation 
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Coen fie denen X 5 (jA1) and Y (after allowance for X,), 

i1.e., find the correlation between the residuals from the 
A 

regression Y = £(X,) and the residuals from another regres- 


1)° 
with Y is now selected and the process continues. As each 


Sion X=f5(X The x. with the highest partial correlation 
variable is entered into the regression, the following 
Values are examined: 

Ee R-, the multiple correlation coefficient; 

2. The partial F-test value for the variable most 
recently entered, which shows whether the variable 
has taken up a significant amount of variation 
over that removed by variables previously entered 
in the regression. 

As soon as the partial F value related to the most recently 


entered variable becomes nonsignificant the process is 


terminated. 
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APPENDIX D 
STATISTICAL TESTS FOR SIGNIFICANCE 


Regression procedures per se may be categorized as de- 
scriptive statistics. Regression analysis is commonly per- 
formed on sample data from which the researcher is either 
interested in estimating population parameters from sample 
regression statistics or to testing statistical hypothesis 
about the population parameters and determining confidence 
limits for estimates in testing the hypothesis. 

The overall test for goodness of fit of the regression 
equation uses statistical inference procedures to test the 
null hypothesis that the sample of observations being analyzed 
has been drawn from a population in which the multiple cor- 
relation 1s equal to zero. An equivalent way of stating 
the null hypothesis, Ho, is that the next variable to be 
added in a forward regression would not add significantly 
to the explained variance in the dependent variable, Y, 
already accounted for by variables included in the regression 
equation. The alternate hypothesis, Hy» directly contradicts 
the null hypothesis. 

The test statistic employed for the overall test is 


SS no 9/K R“/K 


SS_,./(N-K-1) CER oy (Nes 
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where SS reg is the sum of squares explained by the entire 
regression equation, SS... is the residual (unexplained) 

sum of squares, K is the number of independent variables 

in the equation and N is the sample size. The F ratio is 
distributed approximately as the F distribution with degrees 
of freedom K and (N-K-1). 


eZ is an’ R© statistic adjusted for the nunm- 


Adjusted R 
ber of independent variables in the equation and the number 
of cases. It is a more conservative estimate of the percent 
of variance explained, especially when the sample size is 
small. The adjusted R“ formula uses unbiased estimates 


of the error variance and the total variance of Y in the 


pomuiacton, Che formula used by SPSS is: 


Adjusted R® = R* - $=) (1-R*) 


Becaused Adjusted ce is a conservative estimate of the per- 
cent of variance explained, the author used the following 
F statistic to determine confidence levels for ‘the regression 


equations obtained for each command considered: 


(Adjusted R“)/K 
(1-Adjusted oy een 
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